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nighttime  atmosphere,  the  heat  and  humidity  fluxes  exhibited  considerable  variability  in 
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space  and  time.  Consequently,  the  optical  C  ~  calculated  using  flux  measurements  at  a 
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point  showed  considerable  scatter  about  the  path-averaged  optical  measured  by  the 

scintillometer.  Both  incident  short-wave  and  net  radiation  were  found  to  be  convenient 
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indicators  of  temporal  variations  in  optical  CM  during  the  day,  but  not  at  night* 
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Existing  theory  relates  mmw  to  Cj,  ,  C  ,  and  C,j,  .  The  temperature-humidity 


cross-structure  correlation  r„  must  be  known  to  estimate  C„  . 

Tq  Tq 


In  the  absence  of  direct 


measurements  of  r^  ,  it  is  usually  approximated  by  a  constant  between  zero  and  unity. 
Previous  studies  over  moist  soils  have  indicated  that  r,p  is  near  unity,  while  the  present 
DPG  study  over  dry  soil  indicates  that  it  is  near  zero.  A  possible  explanation  for  the  DP< 


results  is  the  lack  of  moisture  available  for  evaporation  from  the  strongly  heated  ground 

2 

surface.  The  previous  tests  of  the  mmw  CM  theory  over  moist  soils  have  confirmed  that  mm 
2  2  N 

C^  is  principally  determined  by  the  contribution,  with  a  secondary  contribution  from 

C^.  However,  in  the  case  of  a  dry  surface,  C.p  can  become  the  dominant  term,  depending  oi 

the  magnitude  of  r^.  Rather  than  make  assumptions  about  the  magnitude  of  ,  it  is 

therefore  preferable  to  obtain  fast-response  (10-Hz)  temperature  and  humidity  measurements 

2 

and  use  spectrum  analysis  techniques  to  determine  r^.  With  r,j,  known,  mmw  can  be 
estimated  using  heat  and  moisture  flux  measurements. 
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SECTION  1 .  SUMMARY 


1 . 1  BACKGROUND 


Interest  in  the  use  of  wave  propagation  techniques  for  military  target 
acquisition  has  been  steadily  increasing  over  the  last  few  decades.  The 
performance  of  sensors  based  on  these  techniques  is  determined  in  part  by  the 
refractive  index  fluctuations  over  the  line  of  sight  between  the  sensor  and 

the  target.  Because  the  effects  of  these  fluctuations  can  be  described  by  the 

2 

refractive  index  structure  parameter  C^  ,  studies  have  been  performed  to 
define  CN^  as  a  function  of  meteorological  variables  such  as  the  temperature 
structure  parameter  C^,  ,  the  humidity  structure  parameter  C^  ,  and  the 
temperature-humidity  structure  parameter  C^.  These  parameters  can  in  turn  be 
related  to  the  atmospheric  fluxes  of  sensible  heat  and  latent  heat  (humidity). 

U.S.  Army  Dugway  Proving  Ground  (DPG)  has  "optical"  scintillometers  that 
2 

are  used  to  measure  CN  for  the  visible  wavelengths,  but  these  measurements 

are  not  made  with  sufficient  regularity  to  develop  an  "optical  climatology" 

for  test  planning  purposes.  Also,  no  instrumentation  is  available  to  measure 
2 

CN  in  other  wavelengths,  in  particular  the  millimeter  wavelengths  (mmw),  of 

special  interest  for  planned  field  tests.  Consequently,  it  is  desirable  to 

2 

establish  working  relationships  between  at  various  wavelengths  and  readily 
measured  variables. 

1 . 2  STUDY  OBJECTIVE 

The  objective  of  the  optical  climatology  study  was  to  define  and  test 

2 

methodologies  for  deriving  optical  and  mmw  C^  from  available  meteorological 
measurements  at  DPG. 

1 . 3  INSTRUMENTATION 

The  principal  instrumentation  used  in  this  study  included  an  optical 
scintillometer,  actinometers,  and  fast-response  vertical  velocity,  tempera¬ 
ture,  and  humidity  sensors.  The  scintillometer  consisted  of  a  transmitter 


1 


with  a  0.94-pm  light  emitting  diode  (LED)  source  and  a  downrange  receiver. 
Scintillation  measurements  were  obtained  as  a  function  of  the  variance  of 
signal  intensity  at  the  receiver.  Single-axis  sonic  anemometer/thermometers 
were  used  to  measure  vertical  velocity  variances  and  the  covariances  of 
temperature  with  vertical  velocity,  trcm  which  sensible  heat  flux  is  derived. 

A  fast-response  hygrometer  was  mounted  alongside  each  sonic  anemometer/ 
thermometer  to  obtain  humidity  variances  and  covariances  of  humidity  with 
temperature  and  vertical  velocity.  Actinometers  included  pyranometers  for 
short-wave  radiation  and  pyrgeometers  for  long-wave  radiation  mounted  to 
measure  incoming  and  outgoing  radiation  at  1.4  m  above  ground  level  (AGL). 

A  2-m  AGL  optical  scintillometer  path  was  selected  at  a  location  repre¬ 
sentative  of  the  DPG  test  grid  terrain.  The  paths  chosen  for  the  summer  and 
winter  trials  differed  because  of  the  requirement  to  keep  direct  sunlight  out 
of  the  scintillometer  optics.  Four  sets  of  the  fast-response  sensors  and  the 
actinometers  were  located  near  the  scintillometer  receiver.  The  fast-response 
instrument  sets  were  placed  in  an  east-west  line  with  6-m  separation  between 
each  set.  Three  of  these  sets  were  mounted  at  3.5  m  AGL  and  one  was  mounted 
at  1.5  m  AGL.  Located  near  this  instrument  cluster  were  two  10-m  meteorologi¬ 
cal  towers  instrumented  for  wind  measurements  at  10  m  and  temperature  and 
humidity  measurements  at  2  m.  Data  were  collected  from  the  fast-response 
instruments  at  10  Hz,  and  from  the  other  instruments  at  1  Hz.  The  1-Hz  data 
were  averaged  into  1-min  data  points  for  storage. 

1.4  RESULTS 

Data  were  collected  during  two  phases:  (a)  Phase  I  in  February  1988  to 

represent  winter  conditions  and  (b)  Phase  II  in  August  1988  to  represent 

summer  conditions.  Many  of  the  desired  ir'-truments  were  unavailable  or 

malfunctioned  during  Phase  I.  Consequent^,  analyses  of  data  from  this  phase 

2 

mainly  addressed  relationships  between  incoming  short-wave  radiation  and  CN  . 
During  Phase  II,  flux  and  variance  data  were  available  from  the  fast-response 
instrumentation  and  were  used  to  derive  C  '  values  for  comparison  with  the 

lx 

values  obtained  from  the  scintillometer  measurements.  Phase  II  actinometer 

2 

data  were  also  compared  with  scintillometer  measurements. 
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2 

The  Phase  II  data  were  used  to  estimate  the  contributions  to  by  the 

fluxes  of  heat  and  humidity  and  the  temperature-humidity  covariance.  Over  a 

strongly  heated  dry  surface  (typical  of  midday  conditions  at  DPG  during  the 

summer),  the  upward  heat  flux  from  the  ground  surface  is  the  major  contributor 
2 

to  optical  CN  .  The  upward  humidity  flux  contributes  only  a  tenth  to  a 
N  2 

hundredth  of  a  percent  to  the  optical  and  can  be  safely  disregarded,  while 

the  temperature-humidity  covariance  term  contributes  between  1  and  5  percent. 

At  night,  the  heat  flux  is  downward  towards  the  surface  while  the  humidity 

flux  remains  upwards.  The  nocturnal  heat  flux  is  10  percent  or  less  of  the 

midday  flux,  and  the  nocturnal  moisture  flux  remains  small,  contributing 
2 

negligibly  to  CN  . 

A  major  difference  between  the  findings  of  this  study  and  previous  stud¬ 
ies  conducted  in  more  humid  conditions  is  that  the  coherence,  or  frequency- 
dependent  correlation  between  temperature  and  humidity,  is  much  less  than 
unity.  Temperature-humidity  coherence  at  midday  tends  to  fall  below  0.75  for 
eddy  sizes  of  10  m  or  less,  due  to  the  lack  of  moisture  available  for  evapor¬ 
ation  at  the  dry  surface.  At  night,  coherence  remains  low  for  all  scales  of 
motion.  Nocturnal  conditions  at  DPG  are  characterized  by  light  winds,  stable 
thermal  stratification,  and  intermittent  turbulence.  The  small  amounts  of 
moisture  available  at  the  surface  often  do  not  rise  through  the  overlying 
layers.  Thermally  stratified  layers  are  subjected  to  differential  rates  of 
cooling,  depending  on  moisture  content  and  contact  with  the  cooler  surface. 
Coherence  is  diminished  as  thermally  isolated  layers,  drifting  about  in  light 
winds,  pass  over  sensors  during  sampling  periods. 

2 

The  optical  CN  computed  for  daytime  conditions  were  consistent  with  and 

"  2  2 

comparable  (within  ±30  percent)  to  measured  CN  .  However,  the  derived 

from  nocturnal  flux  and  covariance  measurements  exhibited  considerable  scatter 

2 

between  instrument  sets  and  differed  from  the  optical  path-averaged 
measurements  by  up  to  several  orders  of  magnitude.  The  intermi t tency  of 

turbulence  in  strongly  stable  nocturnal  conditions  experienced  at  DPG  adds  a 

2 

high  degree  of  uncertainty  to  estimates  of  . 
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usual  assumption  that  the  temperature-humidity  coherence  is  unity,  is  larger 

than  the  temperature  and  humidity  terms.  If  this  assumption  is  correct,  the 

cross-structure  term  nearly  cancels  the  temperature  and  humidity  terms  in  the 

equation,  producing  very  low  mmv  in  dry  desert  conditions.  Analysis 

of  the  DPG  temperature-humidity  cross-spectra  suggests  that  the  temperature- 

humidity  coherence  is  closer  to  zero  than  unity.  These  results  cannot  be 

2 

verified  against  actual  mmv  measurements  because  the  instrumentation 
required  to  obtain  these  measurements  is  not  available  at  DPG. 

2 

Intercomparisons  between  optical  scintillometer  measurements  and 

incoming  short-wave  and  net  radiation  during  both  Phases  I  and  II  indicated 

that  they  covary  during  daylight  hours.  No  relationship  betveen  nocturnal  net 

2 

radiation  and  optical  was  found. 

1 . 5  CONCLUSIONS 

2 

Relationships  between  optical  C^  and  heat  flux  were  verified  for  midday 

over  typical  (summertime)  DPG  test  grid  conditions.  If  a  heat  flux  measure- 

2 

ment  is  available,  it  is  possible  to  calculate  optical  CN  (or  vice  versa) 

2  w 

with  a  high  degree  of  reliability.  Because  CN  can  easily  vary  over  an  order 
of  magnitude  depending  on  the  underlying  ground  conditions,  terrain  varia¬ 
bility  is  the  greatest  source  of  uncertainty  in  the  computation  procedure. 

Incident  short-wave  radiation  or  net  radiation  measurements  can  provide 

2 

valuable  information  on  temporal  variations  in  CN  during  daylight  hours.  The 

2  W 

relationship  between  radiation  and  CN  is  relative  rather  than  absolute 

because  the  magnitude  of  C^  is  also  dependent  upon  ground  conditions  and  wind 

speed.  Nevertheless,  radiation  measurements  provide  a  simple,  inexpensive 

2 

means  to  monitor  CM  and  to  track  it  through  changes  in  cloud  cover. 

N 
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The  relationship  between  path-averaged  measurements  and  point 
measurements  of  heat  flux  at  night  is  tenuous.  The  nocturnal  atm.  phere  is 
highly  stratified,  with  significant  variability  in  space  and  time.  Nocturnal 
turbulence  also  exhibits  a  high  degree  of  intermittency.  For  these  reasons, 
flux  measurements  made  at  a  point  are  not  likely  to  be  representative  of 
conditions  existing  along  an  extended  optical  path.  The  point-to-point 
measurement  variability  at  the  adjacent  fast-response  instrument  sites  was  of 
approximately  the  same  magnitude  as  the  variability  between  the  sites  and  the 
optical  path.  A  more  general  conclusion  is  that  nocturnal  meteorological 
measurements  made  at  one  location  are  not  likely  to  be  representative  of 
conditions  at  other  nearby  locations. 

The  real  component  of  the  mmv  refractive  index  is  composed  of  a 
temperature  structure,  a  humidity  structure,  and  a  temperature-humidity  cross¬ 
structure  term.  The  magnitude  of  the  cross-structure  term  is  dependent  upon 
the  temperature-humidity  coherence.  Previous  mmw  studies  conducted  over  moist 

soils  conclude  that  the  humidity  structure  term  is  the  dominant  contributor  to 
2 

mmw  ,  and  that  temperature-humidity  coherence  is  near  unity.  The  present 
study  indicates  that,  over  dry  soils,  the  temperature-humidity  cross-structure 
term  would  be  dominant  if  it  were  not  diminished  by  a  temperature-humidity 
coherence  much  less  than  unity.  Because  of  the  variability  in  temperature- 
humidity  coherence,  it  is  necessary  to  use  fast-response  temperature  and 
humidity  measurements  and  spectrum  analysis  techniques  to  specify  turbulence 
effects  on  mmw  propagation. 

1 . 6  RECOMMENDATIONS 

The  fast-response  data  sets  collected  during  Phase  II  testing  contain 
roughly  40  megabytes  of  high  quality  data.  Only  a  small  sample  of  the 
available  data  set  was  used  to  generate  the  results  presented  in  this  report. 
Interesting  features  not  addressed  in  this  report  include  the  details  of  flux 
transition  through  sunrise  and  the  presence  of  intermittent  turbulence  bursts 
at  night  that  cause  a  temporary  breakdown  of  the  stability  structure.  Analysis 
of  this  data  set  should  continue. 
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Insufficient  data  were  collected  during  Phase  I  to  analyze  the  relation- 

2 

ship  of  wintertime  fluxes  to  .  The  behavior  of  the  humidity  flux  term  over 
wet  or  snow  covered  surfaces  is  likely  to  be  considerably  different  from  the 
present  dry  ground  results.  These  effects  can  be  particularly  significant  for 
mmw  propagation.  A  second  winter  measurement  project  should  be  initiated 
after  DPG  acquires  its  fast-response  hygrometers  and  a  digital  data  collection 
system. 

2 

Millimeter-wave  cannot  be  calculated  with  certainty  over  dry  soils 
until  the  effects  of  temperature-humidity  coherence  on  the  temperature- 
humidity  cross-structure  term  are  determined.  Co-located  fast-response 
temperature  and  humidity  instruments,  with  alias-free  data  rates  to  10  Hz, 
should  be  used  during  mmw  field  experiments  to  determine  this  coherence,  which 
can  be  applied  as  a  multiplication  factor  on  the  cross-structure  term.  With 
this  modification,  the  mmw  turbulence  theory,  verified  for  moist  conditions 
during  the  Flatville,  Illinois  experiments  (Hill  et  al.,  1988),  can  be  used 
over  dry  soils. 
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SECTION  2.  DETAILS  OF  THE  TEST 


2.1  SITE  DESCRIPTION 

Most  DPG  test  grids  are  located  on  flat  clay  soil  covered  with  a  desert 
shrub  type  vegetation.  The  location  chosen  for  the  optical  climatology  test 
as  being  representative  of  these  test  grids  was  between  the  Atmospheric 
Sciences  Laboratory  (ASL)  rawinsonde  building  (Building  4000)  and  Avery  Area 
(Figure  1).  This  site  has  the  added  convenience  of  readily  available  electric 
power  and  two  10-m  meteorological  masts  instrumented  to  provide  continuous 
information  on  temperature,  relative  humidity,  wind  direction,  and  wind  speed. 
Pressure  and  radiation  data  are  also  continuously  available  at  the  Building 
4000  site. 


Figure  1.  Optical  Climatology  Study  Site. 


The  Building  4000  site  is  located  on  the  north  eastern  corner  of  the 
Ditto  Technical  Center,  several  hundred  meters  south  of  the  eastern  end  of  the 
Michael  Army  Airfield  runway.  There  are  no  significant  topographical  features 
within  1  km  of  the  site  from  the  northwest  quadrant  clockwise  to  the  southeast 
quadrant.  Several  obstructions  (trees,  buildings,  water  tower)  and  a  creek 
bed  are  within  1  km  of  the  site  in  the  southwest  quadrant.  Because  an  unob¬ 
structed  fetch  of  300  to  1000  m  was  desired  for  atmospheric  flux  measurements, 
Phase  II  trials  were  restricted  to  periods  with  either  wind  directions  from 
northeast  to  southeast  or  light  and  variable  winds. 

The  major  difficulty  encountered  at  the  Building  4000  site  was  finding  a 
suitable  scintillometer  path.  The  optical  path  had  to  be  over  relatively 
undisturbed  land  and  fit  within  the  physical  constraints  of  the  runway, 
buildings,  and  roads.  An  additional  constraint  was  that  scintillometer 
transmitter  and  receiver  optics  had  to  be  protected  from  direct  exposure  to 
sunlight,  which  can  occur  when  the  sun  is  at  low  angles  above  the  horizon. 

The  scintillometer  path  was  changed  between  Phases  I  and  II  in  order  to 
protect  the  optics  during  trials  at  sunrise  and  sunset. 

2.2  TEST  DESIGN  AND  INSTRUMENTATION 

The  DPG  instrumentation  used  during  Phase  I  included  a  scintillometer,  a 
three-axis  sonic  anemometer,  the  instrumentation  on  two  10-m  masts,  and  a  set 
of  actinometers  mounted  on  a  platform.  Additional  instrumentation  for  Phase 
II  was  provided  and  operated  by  Campbell  Scientific,  Inc.  (CSI).  These 
sensors  included  four  sets  of  single-axis  sonic  anemometers,  fine-wire 
thermocouples,  and  krypton  hygrometers.  The  three-axis  sonic  anemometer  was 
not  used  during  Phase  II. 

An  optical  scintillometer  (Figure  2),  as  described  in  Ochs  and  Cartwright 

2 

(1980),  was  used  to  measure  the  refractive  index  structure  parameter  C^  .  The 
scintillometer  consists  of  a  transmitter  using  a  0.94-pm  LED  source  and  a  pair 

of  downrange  receivers.  The  log  intensity  variance  in  irradiance  at  the 

2 

receiver  is  used  as  a  measure  of  C^  (see  Appendix  A).  Scintillometer 
measurements  were  taken  at  a  height  of  2  m  above  the  surface.  The  Phase  I 
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path  length  was  674  m,  while  the  Phase  II  path  length  was  195  m.  Power  for 
the  scintillometer  transmitters  was  provided  by  storage  batteries. 


Figure  2.  Model  IV  Optical  Scintillometer  Receiver. 


Actinometer  data  were  taken  from  two  Eppley  precision  spectral 
pyranometers  and  two  Eppley  pyrgeometers ,  as  described  in  the  Eppley 
Laboratory  manual.  These  instruments  were  mounted  on  an  actinometer  platform 
(Figure  3)  to  measure  incident  and  reflected  short-wave  radiation  and  incident 
and  emitted  long-wave  radiation.  The  combination  of  these  measurements 
produced  net  radiation. 
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Figure  3.  DPG  Actinometer  Platform. 


Low  battery  voltages  were  found  for  the  pyrgeometers  at  the  beginning  of 
the  Phase  I  tests.  Thus,  no  reliable  long-wave  radiation  data  were  collected 
until  the  fourth  week.  At  the  end  of  the  Phase  I  test  series,  all  of  the 
actinometers  were  turned  face  up  for  a  quality  control  check.  The  pyrano- 
meters  agreed  to  within  1  percent  of  each  other,  while  the  pyrgeometers  were 
within  3  percent.  The  actinometers  were  sent  to  the  White  Sands  Missile  Range 
calibration  facility  for  calibration  between  Phases  I  and  II.  Unfortunately, 
one  of  the  pyrgeometers  was  damaged  during  return  shipment  and  was  unavailable 
during  Phase  II.  In  order  to  compensate  for  the  missing  pyrgeometer,  incident 
longwave  radiation  was  derived  empirically  using  Brunt's  formula  (Kondratyev, 
1969).  The  procedure  was  similar  to  the  one  used  by  Morgan  et  al.  (1970)  in 
radiation  tests  at  Davis,  California. 
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At  the  conclusion  of  Phase  II  the  remaining  pyrgeometer  was  turned 
face-up  for  2  days  to  allow  comparison  of  the  measured  incident  long-wave 
radiation  with  the  calculated  radiation.  The  linear  correlation  coefficient 
between  measured  and  calculated  radiation  for  Day  1,  which  had  clear  skies 
during  the  entire  day,  was  0.96  and  the  mean  values  agreed  to  within  2 
percent.  Day  2  was  partly  cloudy,  and  the  calculated  radiation  required  an 
adjustment  as  suggested  by  Sellers  (1965).  After  this  adjustment,  the 
correlation  coefficient  was  0.97.  The  means  of  the  calculated  radiation  were 
within  1  percent  of  the  measured  values  for  Day  2.  Because  of  these  favorable 
comparisons,  the  empirical  values  for  incident  long-wave  radiation  were  used 
to  replace  the  missing  pyrgeometer  data  for  net  radiation  computation. 

In  an  attempt  to  obtain  Phase  I  flux  measurements,  a  three-axis  Applied 
Technologies,  Inc.  (ATI)  sonic  anemometer  was  mounted  at  the  2-m  level  on  the 
Base  Weather  Station  10-m  tower.  This  instrument  is  capable  of  measuring  the 
alongwind,  crosswind,  and  vertical  wind  components  and  temperature  for  flux 
measurements.  The  sonic  anemometer  produced  questionable  data  throughout 
Phase  I.  At  times  the  resultant  wind  speed  and  direction  agreed  well  with  the 
tower  data,  while  at  other  times  the  correspondence  was  poor.  The  vertical 
(w)  wind  component  on  the  sonic  anemometer  also  contained  a  persistent  bias, 
suggesting  that  the  instrument  was  not  adequately  leveled.  Consequently,  data 
from  this  instrument  were  not  used  in  the  analysis. 

One  of  the  objectives  of  Phase  II  was  to  evaluate  the  relationships 

2 

between  the  refractive  index  structure  parameter  and  the  measured  sensible 
and  latent  heat  fluxes.  Sensible  and  latent  heat  fluxes  can  be  calculated 
from  the  covariances  of  the  fluctuations  of  the  vertical  wind  speed  (w' ) 
component  with  the  fluctuations  of  the  temperature  (T')  and  absolute  humidity 
(q').  These  variables  were  measured  during  Phase  II  using  single-axis  sonic 
anemometers  with  fine-wire  thermocouples  (Model  CA27T)  and  Krypton  hygrometers 
(Model  KH20)  (see  Tanner  et  al.,  1985).  Four  flux  measurement  systems  were 
set  up,  three  at  a  height  of  3.5  m  AGL  and  one  at  1.5  m  AGL  (Figure  4). 
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Critical  experimental  design  considerations  include  instrument  exposure, 
separation  distance,  path  averaging,  and  data  rate.  The  instruments  were  set 
up  facing  towards  the  southeast,  with  the  hygrometer  probe  set  10  cm  back  from 
and  10  cm  to  the  east  of  the  front  edge  of  the  sonic  anemometer  probe  to  mini¬ 
mize  flow  distortion.  The  fine-wire  thermocouple  probe  was  located  2  to  3  cm 
from  the  sonic  path.  Priestly  and  Hill  (1985)  state  that  sensor  separation 
degrades  the  cospectrum  at  spatial  wavelengths  smaller  than  four  times  the 
sensor  separation.  Consequently,  no  attempt  was  made  to  evaluate  cospectra  at 
scales  smaller  than  50  cm.  Path  averaging  limits  the  spatial  resolution  of 
the  measurement  system.  While  the  spatial  resolution  of  the  thermocouple  and 
hygrometer  is  on  the  order  of  a  few  centimeters,  the  system  resolution  is 
limited  by  the  10-cm  sonic  path  length.  A  data  rate  of  10  Hz  was  found  to  be 
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compatible  with  the  instrument  separation  distance,  resolution,  and  measure¬ 
ment  height. 

2.3  RELATIONSHIP  OF  C  2  TO  OTHER  ATMOSPHERIC  VARIABLES 

N 

2 

Over  the  last  several  decades  a  number  of  studies  have  related  to  the 
fluxes  of  heat  and  humidity.  Coulter  and  Vesely  (1980)  derived  the  empirical 
relationship 


Hs  =  0.48zypcp(g/T)1/2|CNT2/(AP)]3/2  (2-1) 

2 

where  Hg  is  a  derived  heat  flux,  C^  is  the  square  root  of  CN  ,  p  is  air 
density,  cp  is  specific  heat,  T  is  the  mean  temperature,  A  is  a  wavelength- 
dependent  refractive  index  coefficient,  P  is  the  atmospheric  pressure,  z  is 
height  above  the  surface,  and  y  is  a  water  vapor  correction  term.  Equation 
(2-1)  is  strictly  valid  only  for  optical  wavelengths  in  quasi-steady-state, 
convectively  unstable  conditions  such  as  found  at  midday.  At  wavelengths  in 
excess  of  1  um  and  optical  paths  above  warm,  wet  surfaces,  water  vapor 

pressure  fluctuations  become  significant  and  sometimes  emerge  as  the  dominant 

2  2 
contributors  to  CN  .  Vesely  (1976)  consolidated  the  relationship  between  CN 

and  the  fluctuations  of  temperature  and  humidity  into  the  expression 

c„2  ■  <ct2*v/t4,  [i.rTq  m*m2]  <«> 

2 

where  C^,  is  the  temperature  structure  parameter,  B  is  the  Bowen  ratio  (the 
ratio  of  sensible  to  latent  heat  fluxes)  and  r^^  is  the  temperature-humidity 
cross-structure  correlation  coefficient,  given  by 


r 


Tq 


/CTCq 


(2-3) 


Equation  (2-2)  accounts  for  the  combined  temperature  and  humidity  effects  on 
2 

CN  .  These  effects  vary  as  a  function  of  the  covariances  w'T',  w'q',  and 

N  -1  -2  -1 
T'q',  where  w'T'  (K  m  s  )  and  w'q'  (g  m  s  )  determine  respectively  the 

_3 

temperature  and  humidity  fluxes,  and  T'q'  (K  g  m  )  is  the  temperature- 
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humidity  covariance.  Hill  et  al.  (1980)  found  the  pressure  fluctuation  term 
to  be  several  orders  of  magnitude  less  than  the  temperature  and  humidity 
terms.  Consequently,  the  contribution  of  the  pressure  fluctuation  term  is 
ignored  in  Equation  (2-2). 

A  divergence  of  opinion  exists  concerning  the  magnitude  of  r^.  The 
temperature-humidity  cross-structure  correlation  is  dependent  on  the  coherence 
between  T  and  q  at  the  wavelengths  of  turbulence  that  cause  the  scattering 
that  produces  scintillation.  Wesely  and  Hicks  (1978)  and  Coulter  and  Vesely 
(1980)  suggest  that,  at  the  small  spatial  scales  of  refractive  inhomogeneities 
associated  with  the  scattering  of  light,  this  coherence  approaches  zero.  On 
the  other  hand,  Priestley  and  Hill  (1985)  conclude  from  their  measurement 
program  that  T  and  q  are  highly  correlated  at  these  scales,  leading  to  an  r^ 
near  unity.  Gossard  and  Sengupta  (1988)  define  a  homogeneous,  steady-state 
model  in  which  horizontal  gradients  of  variance  and  covariance  can  be 
neglected,  which  allows  r^,^  to  approach  unity  under  ideal  atmospheric  condi¬ 
tions.  The  results  of  the  current  study  indicate  that  the  coherence  between  T 
and  q  is  very  sensitive  to  departures  from  these  ideal  conditions. 

The  turbulence  scale  that  causes  scintillation  response  in  a  scintillome¬ 
ter  is  on  the  order  of  (AA)^"'\  where  A  is  the  scintillometer  operating  wave¬ 
length  and  A  is  the  path  length.  For  an  optical  scintillometer  operating  at  a 
wavelength  of  1  pm  or  less  over  a  1000  m  path,  the  turbulence  scale  is  on  the 
order  of  several  centimeters.  This  turbulence  scale  is  on  the  lower  end  of 
the  inertial  subrange  and  very  near  the  dissipation  range  of  the  spectrum 
where  turbulent  motion  is  dissipated  into  heat.  A  scintillometer  operating  at 
a  wavelength  of  1  cm  over  the  same  path  uses  turbulence  scales  on  the  order  of 
3  m.  Near  the  surface,  this  turbulence  scale  is  likely  to  be  in  the  energy- 
containing  portion  of  the  spectrum  where  turbulence  is  not  isotropic  and 
varying  degrees  of  correlation  exist  between  w,  T,  and  q.  This  subject  is 
further  addressed  in  Subsection  2.4. 

2  2 

Cq  ,  CT  ,  and  C,^  are  related  through  similarity  concepts  to  the  vertical 
fluxes  of  heat  and  humidity.  Free  convection  relations  defined  by  Kohsiek 
(1982)  for  the  unstable  surface  layer  are 
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(2-5) 


v'T'  =  0.55z(g/T)1/2(CT2)3M 

and 

v'q'  =  0.6z(g/T)1/2(CT2)1/4(Cq2)1/2  (2-6) 


2  2 

Kohsiek  (1982)  uses  r_  to  relate  CT  to  C„  and  C 

Tq  Tq  T  q 

third  estimate  of  the  r_,  magnitude  as 

Tq 


Kohsiek  also  offers  a 


_  ...  2A/2.n  2.1/2  ~  n 

rTq  =  CTq/(CT  >  (Cq  >  =  °'75  (2"7) 

Relationships  between  fluxes  and  structure  parameters  are  defined  for 
stable  atmospheric  conditions  using  the  Obukhov  length  L  and  temperature 
scaling  factor  T^.  The  relationship  derived  for  slightly  stable  conditions  by 
Wyngaard  et  al.  (1971)  is 


CT2  =  (T*2/z2/3)[  4.9(1  +  2 . 75(z/L) ) ]  (2-8) 

where  T^  is  the  ratio  of  temperature  flux  w'T'  to  friction  velocity  u^  and  L 
is  defined  by 


L  =  -u*3T/(kgw'T')  (2-9) 

where  k  is  the  Von  Karman  constant.  Nocturnal  conditions  at  DPG  are  charac¬ 
terized  by  clear  skies,  light  winds,  and  extremely  stable  thermal  stratifica¬ 
tion.  These  conditions  violate  assumptions  used  to  develop  T^,  u^,  and  z/L 

relationships.  It  is,  therefore,  unrealistic  to  expect  Equation  (2-8)  to 

2 

define  nocturnal  C^.  with  great  precision. 

Accurate  measurements  of  u^  using  eddy  correlation  techniques  are 
difficult  to  obtain.  Therefore,  the  relationship  between  u*  and  the  squaLe 
root  of  the  variance  in  vertical  wind  speed,  is  used.  Kondo  et  al.  (1978) 
found  that  in  stable  conditions  a  /u.  decreases  from  its  neutral  value  of  1.3, 
reaching  a  plateau  at  1.1  until  the  dimensionless  stability  parameter  z/L 
exceeds  2.0.  For  z/L  greater  than  2.0,  o  /u.  decreases  rapidly  until  it 
stabilizes  near  a  value  of  0.5.  As  a  first  approximation, 


15 


(2-10) 


u .  =  a  / 1.1 

*  w 

offers  values  of  for  the  stability  range  of  interest  in  this  study. 

2.4  DATA  SCREENING  AND  ANALYSES 

Data  from  four  closely  spaced  sets  of  fast-response  (10-Hz)  wind,  temper¬ 
ature,  and  humidity  sensors  were  used  in  the  Phase  II  portion  of  the  study. 

Data  from  these  instruments  were  subjected  to  several  stages  of  screening  prior 
to  analysis.  For  preliminary  screening,  raw  data  were  subjected  to  an  outlier 
check  wherein  data  points  more  than  three  standard  deviations  from  the  mean 
were  flagged.  Very  few  data  points  were  flagged  using  this  procedure,  indi¬ 
cating  that  the  data  did  not  contain  many  large  noise  spikes.  The  data  were 
then  divided  into  18-min  blocks  for  further  reduction  and  analysis  using  the 
FLXVAR  program,  which  calculates  variances  and  covariances,  and  the  RAWSPC 
program,  which  calculates  spectra  and  cross  spectra.  These  programs  permit 
detailed  data  analyses  and  intercomparisons  and  provide  the  variables  used  to 
compute  refractive  index  structure  parameters.  Some  data  sets  were  eliminated 
from  further  analysis  because  of  temperature  drift  in  the  sonic  anemometer 
reference  temperature  junction.  This  screening  procedure  left  16  trials  of 
18-min  duration  for  further  analysis.  Each  data  set  was  further  reduced  by 
selecting  only  data  from  instruments  on  Masts  2  and  3  for  detailed  analysis. 
These  instruments  were  adjacent  (6-m  separation)  and  mounted  at  the  same 
elevation  (3.5  m  AGL). 

Table  1  presents  the  variances  computed  for  vertical  velocity  (w'w'), 
temperature  (T'T'),  and  humidity  (q'q').  The  vertical  velocity  and  tempera¬ 
ture  variances  at  Masts  2  and  3  remain  consistent  with  each  other  for  diurnal 
variations  that  range  over  several  orders  of  magnitude.  The  humidity  variance 
is  uniformly  small  due  to  the  dry  ground  conditions  10  days  after  the  last 
major  rainfall.  Comparison  of  the  humidity  variances  for  the  hygrometers  on 

Masts  2  and  3  also  reveals  a  measurement  bias,  with  the  magnitudes  of  the  Mast 

2  -6 

3  humidity  variances  approximately  0.0175  g  m  greater  than  those  of  Mast  2. 
Analysis  of  the  spectra  for  data  from  these  instruments  revealed  a  large  noise 
component  on  Mast  3  at  the  high  frequency  end  of  the  measured  spectra,  and  the 
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Mast  3  humidity  variance  bias  is  thought  to  be  due  to  this  noise.  The  noise 
effects  did  not  appear  in  covariances  because  the  random  noise  signal  is 
uncorrelated  with  temperature  and  vertical  velocity. 


Table 

1.  Variances  of 

Vertical 

Veloci ty 

(w'w' ), 

Temperature 

(T'T' ) , 

and 

Humidity  (q'q'). 

Start 

Mast  2 

Mas  t  3 

Mast  2 

Mast  3 

Mast  2 

Mast  3 

Trial 

Date 

Time 

w' w' 

w'  w' 

•p  t  'p  t 

•p  t  *p  t 

q'q' 

q'q' 

(1988) 

(MDT) 

✓  2  -2, 
(m  s  ) 

✓  2  ~2\ 
(m  s  ) 

2 

(K  ) 

2 

(K  ) 

,  2  -6, 
(g  m  ) 

.  2  -6v 

(g  m  ) 

1A 

16  Aug 

1105 

.121 

.133 

.699 

.525 

.014 

.027 

IB 

16  Aug 

1124 

.119 

.144 

.534 

.610 

.009 

.027 

ID 

16  Aug 

1159 

.136 

.15  ■ 

.640 

.642 

.010 

.029 

2A 

17  Aug 

0356 

.002 

.002 

.055 

.059 

.015 

.033 

2B 

17  Aug 

0415 

.012 

.012 

.085 

.074 

.052 

.062 

2C 

17  Aug 

0435 

.012 

.010 

.065 

.070 

.002 

.022 

3A 

17  Aug 

0621 

.003 

.003 

.165 

.166 

.014 

.033 

3B 

17  Aug 

0640 

.010 

.010 

.731 

.762 

.043 

.063 

AC 

17  Aug 

1240 

.146 

.142 

.612 

.580 

.082 

.098 

4D 

17  Aug 

1259 

.185 

.186 

.794 

.759 

.017 

.033 

5A 

18  Aug 

0421 

.015 

.012 

.510 

.504 

.009 

.025 

5B 

18  Aug 

0440 

.006 

.006 

.187 

.172 

.006 

.025 

6A 

18  Aug 

0631 

.002 

.002 

.167 

.158 

.015 

.034 

6B 

18  Aug 

0650 

.007 

.006 

.800 

.693 

.021 

.035 

6C 

18  Aug 

0710 

.006 

.005 

.163 

.156 

.009 

.028 

7A 

18  Aug 

1101 

.104 

.121 

.464 

.566 

.013 

.036 

The  temperature  and  vertical  velocity  covariances  at  Masts  2  and  3 
exhibited  large  positive  values  at  midday  and  small  negative  values  in  the 
early  morning  hours.  These  results,  which  are  shown  in  Table  2,  are 
consistent  in  sign  and  magnitude  with  the  results  obtained  in  the  1968  Kansas 
Field  Program  (Izumi,  1971). 
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The  covariances  of  humidity  with  vertical  velocity  were  uniformly  posi¬ 
tive,  indicating  continued  evaporation  from  the  ground  surface  following  10 
days  of  intense  heating  after  the  last  recorded  precipitation.  The  covari¬ 
ances  measured  at  night  and  during  the  transition  at  sunrise  were  character¬ 
ized  by  occasional  bursts  of  flux  interspersed  with  quiescent  periods  of  very 
low  turbulence.  Covariances  during  these  periods,  which  were  two  or  three 
orders  of  magnitude  below  midday  covariances,  were  so  small  that  they  could 
have  been  influenced  by  neglected  secondary  effects  such  as  advection.  Con¬ 
sequently,  they  were  assigned  a  magnitude  of  0.000  in  Table  2. 

The  data  in  Table  2  were  used  to  compute  the  Bowen  ratio,  the  ratio  of 
sensible  to  latent  heat  flux  densities,  using  the  expression 

c  P(w'T')  (2-11) 

B  =  -2— - 

eL  (w'e') 
w 

where  t  is  the  ratio  of  the  molecular  weights  of  water  vapor  to  dry  air,  Ly  is 
the  latent  heat  of  vaporization,  and  w'e'  is  the  humidity  flux  in  units  of 
millibar  meters  per  second  (mb  m  s Vapor  pressure  e  in  millibars  is 
related  to  absolute  humidity  q  in  grams  per  cubic  meter  and  air  temperature  T 
in  degrees  Kelvin  by 


qT  (2-12) 

e  =  - 

216.7 

Table  3  summarizes  the  Bowen  ratios  calculated  using  Equation  (2-11). 

The  Bowen  ratios  B  for  the  daytime  trials  range  between  1  and  10,  indicating  a 
strong  positive  (upward)  heat  flux  and  a  weak,  upward  humidity  flux  above 
strongly  heated  dry  ground.  At  night,  B  changes  sign  due  to  the  downward  heat 
flux.  The  smaller  heat  flux  at  night  occasionally  caused  the  absolute  value 
of  B  to  be  slightly  less  than  1.0.  For  those  cases  shown  in  Table  3  in  which 
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Table  2.  Covariances  of  Vertical  Velocity  and  Temperature  (v'T'),  Vertical 
Velocity  and  Humidity  (w'q'),  and  Temperature  and  Humidity  (T'q'). 


Start 

Mast  2 

Mast  3 

Mast  2 

Mast  3 

Mast  2 

Mast  3 

Trial 

Date 

Time 

w'T' 

w'T' 

w'q' 

w'q' 

T'q' 

T'q' 

(1988) 

MDT 

(K  m  s1) 

(K  m  s-1)(g 

-2  -1 
m  s 

)(g  m"2  s" 

S(K  g  m-3 

)(K  g  nf3) 

1A 

16  Aug 

1105 

.156 

.115 

.018 

.015 

.087 

.066 

IB 

16  Aug 

1124 

.134 

.150 

.013 

.016 

.057 

.069 

ID 

16  Aug 

1159 

.159 

.164 

.016 

.018 

.057 

.064 

2A 

17  Aug 

0356 

-.001 

.000 

.000 

.000 

-.006 

-.006 

2B 

17  Aug 

0415 

-.006 

-.006 

.003 

.002 

-.040 

-.035 

2C 

17  Aug 

0435 

-.005 

-.005 

.000 

.000 

-.002 

-.002 

3A 

17  Aug 

0621 

-.004 

-.003 

.001 

.002 

-.030 

-.033 

3B 

17  Aug 

0640 

-.010 

-.014 

.003 

.004 

-.149 

-.156 

4C 

17  Aug 

1240 

.139 

.120 

.011 

.013 

.062 

.063 

4D 

17  Aug 

1259 

.176 

.156 

.025 

.024 

.086 

.091 

5A 

18  Aug 

0421 

-.030 

-.015 

.004 

.003 

-.048 

-.055 

5B 

18  Aug 

0440 

-.007 

-.006 

.002 

.000 

-.017 

-.015 

6A 

18  Aug 

0631 

-.001 

.001 

.000 

.000 

-.044 

-.044 

6B 

18  Aug 

0650 

-.012 

-.002 

.002 

.000 

-.115 

-.106 

6C 

18  Aug 

0710 

-.003 

.000 

.000 

.000 

-.016 

-.018 

7A 

18  Aug 

1101 

.110 

.116 

.016 

.027 

.056 

.083 

the  humidity  flux  and  B  are  missing,  the  humidity  flux  was  below  measurement 
resolution.  It  is  reasonable  to  assume  in  these  cases  that  B  is  a  large 
number  of  undetermined  sign.  As  a  first-order  approximation,  the  effects  of 
humidity  as  represented  by  B  in  the  denominator  of  the  moisture  contribution 
terms  of  Equation  (2-2)  can  oe  neglected  over  strongly  heated  dry  ground 
during  summer  months  at  DPG. 


Before  calculating  ,  it  is  necessary  to  define  the  energy  containing 
frequencies  in  the  data  and  determine  the  magnitude  of  the  temperature- 
humidity  cross-structure  correlation  coefficient  r,p  .  This  parameter  is  a 
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Table  3.  Mean  Wind  Speeds  (WS),  Temperatures  (T),  Pressures  (P),  Latent  Heats 
of  Vaporization  (L  ),  Humidity  Fluxes  (w'e'),  and  Bowen  Ratios. 

W 


Trial 

Date 

(1988) 

Start 

Time  WS 
(MDT) (m  s-1 

T 

)(°C) 

P 

(mb) 

L 

-1 

(c  g  ) 

Mast  2 

v'e' 

(mb  m  s  1) 

Mast  3 

v'e' 

(mb  m  s  1) 

Bowen 

Mast 

2 

Ratio 

Mast 

3 

1A 

16 

Aug 

1105 

3.0 

25.4 

869.5 

583 

.026 

.021 

3.5 

3.2 

IB 

16 

Aug 

1124 

3.0 

26.0 

869.5 

583 

.018 

.023 

4.3 

3.7 

ID 

16 

Aug 

1159 

2.6 

27.0 

869.4 

582 

.022 

.025 

4.2 

3.8 

2A 

17 

Aug 

0356 

2.2 

12.8 

867.8 

590 

a 

— 

— 

-- 

2B 

17 

Aug 

0415 

2.2 

11.8 

867.8 

590 

.004 

.003 

-0.9 

-1.2 

2C 

17 

Aug 

0435 

1.8 

11.1 

867.9 

591 

— 

— 

— 

— 

3A 

17 

Aug 

0621 

1.8 

10.1 

868.4 

592 

.002 

.003 

-1.2 

-0.7 

3B 

17 

Aug 

0640 

1.8 

10.3 

868.4 

592 

.004 

.005 

-1.4 

-1.5 

4C 

17 

Aug 

1240 

1.7 

33.8 

868.6 

578 

.016 

.019 

5.0 

3.7 

4D 

17 

Aug 

1259 

3.2 

33.7 

868.4 

578 

.036 

.034 

2.9 

2.7 

5A 

18 

Aug 

0421 

2.8 

13.6 

867.1 

589 

.005 

.004 

-3.7 

-2.3 

5B 

18 

Aug 

0440 

2.6 

12.9 

867.1 

589 

.002 

— 

-1.9 

-- 

6A 

18 

Aug 

0631 

1.4 

13.2 

867.8 

589 

— 

— 

— 

— 

6B 

18 

Aug 

0650 

1.6 

13.0 

867.8 

589 

.002 

-- 

-3.0 

— 

6C 

18 

Aug 

0710 

1.6 

12.9 

867.9 

589 

— 

— 

— 

— 

7  A 

18 

Aug 

1101 

1.4 

34.1 

869.0 

578 

.023 

.038 

2.8 

1.8 

g 

Humidity  flux  insignificantly  different  from  zero. 


function  of  scintillation  wavelength,  and  is  therefore  a  function  of  fre¬ 
quency.  Spectrum  analysis  provides  a  breakdown  by  frequency  band  of  contri¬ 
butions  to  the  total  variance  in  the  data.  Spectra  of  the  vertical  velocity, 
temperature,  and  humidity  data  for  each  trial  were  examined  using  the  DPG 
spectrum  analysis  program  described  in  Appendix  B.  The  10-Hz  temperature  and 
humidity  data  points  were  instantaneous  readings  and  contained  no  averaging 
from  higher  frequency  data.  Therefore,  the  spectra  contained  energy  aliased 
from  the  unsampled  higher  frequencies.  To  minimize  these  aliasing  effects, 
1-hr  blocks  of  10-Hz  data  were  averaged  over  a  frequency  decade  to  1-s  inter¬ 
vals  prior  to  analysis  of  the  spectra.  This  averaging  virtually  eliminated 
aliasing  effects  that  were  evident  in  unaveraged  10-Hz  data. 
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Spectral  frequencies  were  converted  to  wavelength  by  dividing  the  trial- 
averaged  3.5-m  windspeed  by  the  central  frequency  of  each  spectral  frequency 
band.  Wind-speed  measurements  at  10  m  were  adjusted  to  3.5  m  using  simple 
power-law  relationships  defined  by  Biltoft  (1978).  The  wavelengths  obtained 
for  each  of  the  1-hr  blocks  of  1-Hz  data  ranged  from  several  meters  at  the 
high  frequency  end  of  the  spectrum  to  nearly  1000  m  at  the  low  frequency  end. 

Review  of  the  spectra  also  served  as  an  important  quality  control  check 
to  identify  high  frequency  noise  contamination  independent  of  the  aliasing 
problem.  Unlike  aliasing,  the  noise  was  not  eliminated  by  averaging.  High 
frequency  noise  was  identified  in  the  Mast  3  humidity  data,  which  caused 
higher  variances  in  this  unit  than  those  obtained  from  the  other  humidity 
instruments.  Because  this  noise  was  uncorrelated  with  the  temperature  and 
vertical  velocity  data,  the  covariances  presented  in  Table  2  remain 
unaffected. 

Spectra  of  surface  layer  micrometeorological-scale  atmospheric  variables 
can  be  divided  roughly  into  three  regions  (Kaimal,  1988):  (a)  the  energy 

containing  range,  where  energy  is  produced  by  buoyancy  and  shear;  (b)  the 
inertial  subrange  where  energy  is  neither  produced  nor  dissipated,  but  cas¬ 
cades  progressively  towards  smaller  scales;  and  (c)  the  dissipation  range 
where  kinetic  energy  is  converted  to  internal  energy  through  viscous  dissipa¬ 
tion.  The  portions  of  the  spectra  in  the  energy  containing  range  include 
scales  ranging  in  wavelength  from  meters  to  hundreds  of  meters.  Covariances  of 
variables  in  the  energy  containing  range  are  primarily  responsible  for  the 
fluxes  of  heat,  momentum,  and  water  vapor. 

The  inertial  subrange  includes  scales  ranging  from  a  few  millimeters  to 
several  meters.  Within  the  convective  surface  boundary  layer,  the  outer 
scale,  or  low  frequency  end  of  the  inertial  subrange,  is  approximately  equal 
to  the  height  above  the  surface.  Energy  within  the  inertial  subrange  is 
transferred  in  an  essentially  loss-free  manner  to  progressively  smaller  scales 
and  becomes  partitioned  more  or  less  equally  in  all  spatial  dimensions.  As 
this  occurs,  covariances  between  vertical  motion  and  transported  atmospheric 
constituents  such  as  heat  and  water  vapor  diminish.  Energy  cascading  down  the 
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scales  in  the  inertial  subrange  finally  reaches  the  dissipation  range,  where 
viscous  effects  become  dominant.  At  scales  of  a  millimeter  and  below,  random 
turbulent  motions  are  converted  into  molecular  motions  or  heat. 

The  maximum  covariance  between  w,  T,  and  q  can  occur  when  they  vary  over 
the  same  range  of  scale.  Consequently ,  the  vertical  velocity,  temperature,  and 
humidity  spectra  were  examined  to  determine  the  range  of  wavelengths  for  spec¬ 
tral  maxima  as  a  function  of  time  of  day.  The  ranges  of  spectral  maxima  for 
midday  and  nocturnal  conditions  ate  presented  in  Table  4.  For  trials  conducted 
at  midday,  Table  4  shows  that  the  spectral  maxima  for  w,  T,  and  q  are  of  simi¬ 
lar  magnitude  and  are  well  within  the  energy  containing  region  of  the  spectrum. 
Consequently,  a  relatively  high  degree  of  covariance  was  found,  which  results 
in  large  average  flux  values  (see  Table  5).  At  night,  the  wavelength  scales  of 
the  spectral  maxima  for  w  are  significantly  smaller  than  the  midday  scales, 
while  the  spectral  maxima  for  T  and  q  occur  at  longer  wavelengths.  This  de¬ 
crease  in  the  vertical  motion  scale  is  due  less  to  a  shift  in  spectral  power 
peaks  than  to  a  differentially  greater  loss  of  power  at  the  longer  wavelengths. 
The  shift  of  the  w  maximum  towaid  the  lower  end  of  the  energy  containing  range 
of  the  spectrum,  and  the  attendant  scale  separation  of  spectral  maxima  is 
largely  responsible  for  the  smaller  fluxes  of  heat  and  humidity  during  the 
night . 

Table  4.  Wavelength  Ranges  of  Spectral  Maxima  For  Vertical  Velocity, 
Temperature,  and  Humidity. 

Time  Wavelength  Range  of  Spectral  Maximum  (m) 

of  - 

Day  Vertical  Velocity  w  Temperature  T  Humidity  q 


Midday 

10  -  300 

20  - 

300 

30  - 

300 

Night 

5  -  40 

300  - 

1000 

250  - 

1000 

The  averaged  total  heat  and  humidity  fluxes  and  the  flux  contributions 
from  the  region  of  spectral  maxima  are  shown  in  Table  5.  Comparison  of  the 
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midday  flux  contribution  from  the  spectral  maximum  with  the  total  flux  illus¬ 
trates  the  large  contribution  from  a  rather  narrow  band  within  the  energy 
containing  region  of  the  spectrum.  At  night  and  during  the  sunrise  transi¬ 
tion,  the  contributions  to  the  total  fluxes  are  smaller  and  more  distributed. 


Table  5.  Total  And  Spectral  Maxima  For  Heat  (w'T')  And  Humidity  (w'q')  Fluxes. 

Total  (V  m  2) 

Spectral  Maxima  (V 

m  2) 

Day 

Night 

Sunrise 

Day 

Night 

Sunrise 

w'T'  w'q' 

w'T'  w'q' 

w'T'  w'q' 

w'T'  w'q' 

w'T'  w'q' 

w'T'  w'q' 

142  45 

-8  3 

-5  3 

48  15 

-2  0.8 

-0.5  0.7 

The  w'T'  and  w'q'  fluxes  indicate  the  movement  of  heat  and  humidity  verti¬ 
cally  through  the  atmosphere.  The  larger  scale  movements  or  fluxes  generated 
within  the  energy  containing  range  of  the  spectrum  transport  smaller  scale 
scintillation-causing  eddies  across  an  optical  path.  Therefore,  the  flux 

summed  across  the  entire  energy  containing  range  is  the  pertinent  variable  for 

2 

w'T'  and  w'q'  effects  on  .  In  contrast,  the  T'q'  covariance  is  pertinent 
only  for  the  wavelength  bands  that  affect  the  scintillation  scales  of  interest. 
As  discussed  in  Subsection  2.3,  the  wavelength  bands  of  interest  range  from 
several  millimeters  to  several  centimeters  for  optical  wavelength  scintillation 
effects,  and  over  several  meters  for  mmw  scintillation  effects.  Figures  5  and 
6  show  the  wavelength-dependent  correlation,  or  coherence,  between  combinations 
of  w' ,  T' ,  and  q'  for  summer  midday  (Trial  1)  and  nocturnal  (Trial  2)  condi¬ 
tions  at  DPG.  Because  of  the  averaging  done  to  eliminate  aliasing,  coherence 
in  Figures  5  and  6  does  not  extend  to  wavelengths  much  below  10  m.  However,  a 
trend  towards  decreasing  coherence  at  shorter  wavelengths  is  apparent  in  these 
Figures. 

The  w'T'  and  w'q'  coherence  is  relatively  high  within  the  energy  contain¬ 
ing  region  of  the  spectrum  during  midday.  This  coherence  is  low  at  night  due 
largely  to  the  weakness  of  vertical  motions  within  the  energy  containing 
region  of  the  spectrum.  There  is  also  a  trend  for  coherence  to  decrease  at 
progressively  smaller  wavelengths  within  the  inertial  subrange  due  to  random 
partitioning  of  turbulent  energy. 
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Coherence  of  v'  to  q',  T'  to  q'  and  w'  to  T'  for  Trial  1,  Repre¬ 
senting  Midday  Conditions  During  the  Summer  at  DPG. 
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Figure  6.  Coherence  of  v'  to  q',  T'  to  q',  and  w'  to  T'  for  Trial  2,  Repre 
senting  Nocturnal  Conditions  During  the  Summer  at  DPG. 
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The  midday  T'q'  coherence  shown  in  Figure  5  remains  high  for  wavelengths 
in  excess  of  10  m,  but  decreases  at  shorter  wavelengths.  The  nocturnal  T'q' 
coherence  in  Figure  6  remains  low  for  all  frequencies.  Although  the  spectra 
do  not  extend  to  the  centimeter  scales  to  which  optical  scintillometers 
respond,  trends  in  the  available  data  suggest  that  coherence  at  these  scales 
is  closer  to  zero  than  to  unity.  These  T'q'  coherence  tendencies  differ  from 
the  findings  of  Priestley  and  Hill  (1985),  who  reported  coherence  of  0.8  or 
greater  extending  to  wavelengths  well  below  1  m  for  measurements  over  moist 
soils . 

Several  possible  explanations  exist  for  the  low  T'q'  coherence  in  the  DPG 
data.  One  possibility  is  an  undetected  deficiency  in  data  collection  or 
processing.  However,  the  data  reduction  computer  program  used  on  the  DPG  data 
was  derived  from  the  one  described  by  Kaimal  and  Gaynor  (1983)  and  should  be 
appropriate  for  use  with  these  data.  Also,  the  10-Hz  data  were  averaged  to  1 
Hz  to  eliminate  aliasing,  and  these  averaged  data  exhibited  no  obvious 
remnants  of  aliasing  that  would  artificially  degrade  coherence. 

Another  possible  explanation  for  the  low  T'q'  coherence  in  the  DPG  data 
lies  within  the  physics  of  the  problem.  Neglecting  molecular  diffusion  and 
assuming  that  a  flux  tends  to  be  directed  down  the  gradient  (from  the  source 
towards  a  sink),  high  T'q'  coherence  can  occur  if  the  ratio  of  moisture  and 
heat  gradients  remains  constant.  Gossard  and  Sengupta  (1988)  represent  this 
ratio  as  ( 3q/3z)/( 30/3z) ,  where  0  is  potential  temperature.  They  conclude 
that  this  ratio  is  a  conservative  property  over  homogeneous  terrain,  causing 
high  T'q'  coherence  down  to  very  small  spatial  scales.  While  the  ratio  of 
gradients  for  saturation  vapor  pressure  to  temperature  ( 3es/3z)/( 3T/3z)  is 
conservative,  the  actual  vapor  pressure  or  absolute  humidity  above  the  ground 
surface  is  influenced  by  factors  such  as  percolation  of  moisture  through  the 
soil  and  variations  in  evapotranspiration  rates.  As  long  as  sufficient 
moisture  is  available  at  the  surface,  the  conclusion  of  Gossard  and  Sengupta 
(1988)  remains  approximately  correct.  However,  as  a  surface  becomes 
progressively  drier,  less  moisture  is  available,  and  the  parts  of  the  surface 
that  have  been  most  strongly  heated  are  likely  to  have  the  least  available 
moisture.  The  result  is  a  decrease  in  T'q'  coherence.  Coherence  should 
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approach  zero  when  the  moisture  gradient  vanishes  above  a  desiccated  surface. 
The  lowest  coherence  (0.76)  reported  by  Priestly  and  Hill  (1985)  occurred 
during  their  driest  daytime  test  conditions,  a  result  that  is  consistent  with 
this  analysis. 


An  additional  complication  occurs  during  measurements  at  night  when  a 
stable  thermal  stratification  isolates  the  surface  layer  from  layers  several 
meters  above  it,  eliminating  convective  transfer  of  heat  and  moisture.  These 
upper  layers  of  air  are  drifting  randomly  over  the  surface,  mixing  very  little 
with  layers  above  and  below.  Also,  shifts  in  wind  direction  may  cause  instru¬ 
ments  at  a  fixed  location  to  sample  air  from  several  different  origins  over 
the  course  of  a  measurement  period.  Therefore,  the  T'q'  coherence  for  air 
sampled  during  dry  nocturnal  conditions  is  likely  to  have  a  large  random  com¬ 
ponent  that  differs  significantly  from  unity,  as  demonstrated  in  Figure  6. 


2.5  CALCULATION  OF  OPTICAL  AND  MILLIMETER-WAVE  CM  FROM  FLUXES 

N 

2 

Equation  (2-2)  can  be  used  to  compute  optical  CN  from  the  fluxes  of  heat 

and  moisture  as  measured  by  the  fast-response  instruments  on  Masts  2  and  3. 

Intermediate  steps  in  this  procedure,  which  include  the  estimation  of  the 

2 

Bowen  ratio  and  r_  ,  are  discussed  in  Subsection  2.4.  Estimates  of  C„  and 
2  Tq  T 

Cq  are  obtained  in  this  subsection  using  the  method  of  Kohsiek  (1982)  for  the 
unstable  trials  and  the  method  of  Wyngaard  et  al.  (1971)  for  the  stable 
trials.  Because  the  trials  conducted  during  the  sunrise  transition  period 
generally  had  negative  average  heat  fluxes,  the  method  of  Wyngaard  (1971)  was 
used  for  these  trials  as  well. 


2  2 

The  C,j,  and  C^  values  calculated  using  Equations  (2-5)  and  (2-6)  for  the 

convectively  unstable  midday  data  found  in  Tables  2  and  3  are  shown  in  Table 

6.  These  data  were  computed  for  the  2-m  level,  assuming  constant  fluxes 

between  2  m  and  the  3.5-m  measurement  level.  The  relative  magnitudes  of  the 
2  2 

CT  and  C  values  in  Table  6  indicate  the  relative  strengths  of  daytime  heat 
and  humidity  fluxes  over  dry  desert  terrain.  When  presented  in  terms  of  Bowen 
ratio,  the  CT  contribution  (represented  by  the  term  in  Equation  (2-2)  con- 
taining  r^  )  contributes  between  1  and  3  percent  to  the  total  optical  C^ 


26 


if  an  r_  of  1.0  is  assumed. 
Tq 


Because  r,^  at  optical  wavelengths  is  likely  to 

be  closer  to  zero,  this  contribution  becomes  negligible.  The  C  ^  contribu- 

2  ^ 

tion,  represented  by  the  (.0313/B)  term  in  Equation  (2-2),  is  on  the  order  of 


0.01  percent,  which  is  also  negligible. 


2  2 

Table  6.  Calculated  Temperature  (C„  )  and  Humidity  (C  )  Structure  Parameters 
for  the  Midday  Trials.  ^ 


Trial 


Date 

(1988) 


Start 

Time 

(MDT) 


,2  , 


2/3 


Un 


v  -3,2., 

(g  m  ) 


2/3 


Mast  2  Mast  3 


Mast  2  Mast  3 


1A 

16  Aug 

1105 

.724 

.480 

.009 

.007 

IB 

16  Aug 

1124 

.590 

.683 

.005 

.007 

ID 

16  Aug 

1159 

.742 

.772 

.006 

.008 

4C 

17  Aug 

1240 

.632 

.519 

.004 

.005 

AD 

17  Aug 

1259 

.863 

.734 

.015 

.015 

7A 

18  Aug 

1101 

.463 

.499 

.009 

.022 

2  2 

The  C„  and  C  values  calculated  for  the  nocturnal  and  stable  transition 
T  q 

period  trials  using  Equation  (2-8)  are  given  in  Table  7.  Following  Hill 
2 

et  al.  (1988),  C  was  calculated  by  replacing  T.  in  Equation  (2-8)  with  the 

q 

humidity  scaling  factor  q^,  which  is  the  ratio  of  w'q'  to  u^.  The  use  of 

Equation  (2-8)  required  as  initial  steps  the  estimation  of  u^  from  Equation 

(2-10)  and  the  Obukhov  length  L  from  Equation  (2-9).  As  in  the  case  of  Table 
2 

6,  C,j.  was  computed  for  the  2-m  level  by  assuming  a  constant  flux  between  that 

level  and  the  instrument  height  of  3.5  m.  This  assumption  is  not  well  founded 

2 

for  extremely  stable  conditions.  The  nocturnal  C^.  in  Table  7  are  generally 

smaller  than  the  daytime  values  in  Table  6.  Also,  the  differences  between  the 
2 

CT  calculated  using  data  from  Masts  2  and  3  are  sometimes  large,  indicating 
the  general  lack  of  homogeneity  in  very  stable  nocturnal  conditions.  No 
computations  were  made  using  Mast  3  data  for  Trials  2A,  6A,  and  6C  because  the 
signs  of  the  heat  fluxes  for  these  periods  were  indeterminate. 


Table  7.  Friction  Velocity  (u*),  Obukhov  Length  (L),  and  Calculated 

2  2 

Temperature  (C^,  )  and  Humidity  (C  )  Structure  Parameters  for  the 
Stable  Trials. 


Trial 

Start 

u*(m 

s"1) 

L(m) 

C  2  - 
°T 

b 

K2 

2/3 

n 

r/  -3.2, 

,  (g  m  ) 

c  c 

q  2/3 

m 

(1988) 

(MDT) 

Mast  2 

Mast  3 

Mast  2 

Mast  3 

Mast  2 

Mast  3 

Mast  2 

Mast  3 

2A 

17  Aug 

0356 

.04 

Ma 

7.2 

M3 

.002 

Ma 

Ma 

Ma 

2B 

17  Aug 

0415 

.10 

.10 

11.5 

12.3 

.018 

.016 

.005 

.002 

2C 

17  Aug 

0435 

.10 

.09 

14.5 

14.5 

.011 

.012 

Ma 

Ma 

3A 

17  Aug 

0621 

.05 

.05 

2.8 

2.6 

.049 

.039 

.006 

.015 

3B 

17  Aug 

0640 

.09 

.09 

5.8 

3.7 

.065 

.196 

.006 

.016 

5A 

18  Aug 

0421 

.11 

.10 

3.4 

4.9 

.585 

.149 

.008 

.005 

5B 

18  Aug 

0440 

.07 

.07 

3.4 

4.1 

.087 

.050 

.005 

Ma 

6A 

18  Aug 

0631 

.04 

Ma 

7.2 

Ma 

.002 

Ma 

Ma 

Ma 

6B 

18  Aug 

0650 

.07 

.07 

2.5 

11.3 

.266 

.005 

.005 

Ma 

6C 

18  Aug 

0650 

.07 

Ma 

8.1 

Ma 

.008 

Ma 

Ma 

Ma 

a  Data  Unavailable. 

2 

Table  8  compares  the  values  calculated  using  Equation  (2-2)  with  the 

2 

values  measured  by  the  optical  scintillometer.  The  calculated  C„  values  in 

2  2  4  N 

this  table  vary  according  to  the  multiplication  factor  (A'r/T4 ),  which  repre¬ 
sents  the  sensitivity  of  the  refractive  index  for  the  optical  wavelengths  to 

the  average  pressure  and  temperature.  The  multiplication  factors  are  included 

2 

in  Table  8.  Most  of  the  values  computed  for  the  midday  trials  are  within 

±30  percent  of  the  observed  values.  This  degree  of  correspondence  is  indica- 

tive  of  a  well  mixed  atmosphere.  Although  the  CN  values  for  some  of  the 

nocturnal  and  transition  trials  agree  well  with  the  observed  values,  others  are 

considerably  different.  This  hit-or-miss  condition  is  likely  due  to  localized 
2 

differences  in  CN  that  occur  during  nocturnal  stratification  and  intermittent 
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turbulence  conditions.  The  overall  precision,  defined  as  the  standard  devia- 

2 

tion  of  the  differences  between  observed  and  measured  C„  values  for  all  tri- 

N 

als,  is  about  an  order  of  magnitude.  This  result  is  strongly  influenced  by  one 

2  2 
calculated  (Trial  5A,  Mast  2)  that  diverged  widely  from  the  measured  CN  . 

For  the  nocturnal  and  transition  trials,  the  precision  of  the  C^2  calculated 

using  the  Mast  2  and  Mast  3  data  was  roughly  comparable  to  the  precision  ob- 

2 

tained  when  the  CN  calculated  using  data  from  either  mast  are  compared  to  the 

2  N  2 
measured  Cjj  ,  further  illustrating  the  effects  of  localized  variations  on  . 


Table  8.  Multiplication  Factor  and  Calculated  and  Measured  Optical  Refractive 
Index  Structure  Parameters.3 


Trial 

Date 

(1988) 

Start 

Time 

(MDT) 

Multiplication 

Factor 

?  2  4 

(AP/r) 

2 

Calculated  C„ 

N 

Mast  2  Mast  3 

Measured 
C  2 

N 

1A 

16  Aug 

1105 

57.76 

41.8 

27.7 

50.6 

IB 

16  Aug 

1124 

57.76 

34.1 

39.5 

56.2 

ID 

16  Aug 

1159 

56.25 

41.7 

43.4 

59.4 

2A 

17  Aug 

0356 

68.89 

0.1 

Mb 

4.1 

2B 

17  Aug 

0415 

70.56 

1.3 

1.1 

2.2 

2C 

17  Aug 

0435 

70.56 

0.8 

0.8 

1.9 

3A 

17  Aug 

0621 

72.25 

3.5 

2.8 

5.7 

3B 

17  Aug 

0640 

72.25 

4.7 

14.2 

6.0 

4C 

17  Aug 

1240 

51.84 

32.8 

26.9 

34.7 

4D 

17  Aug 

1259 

51.84 

44.7 

38.1 

41.2 

5A 

18  Aug 

0421 

67.24 

39.3 

10.0 

6.4 

5B 

18  Aug 

0440 

68.89 

6.0 

3.4 

4.3 

6A 

18  Aug 

0631 

68.89 

0.1 

Mb 

2.6 

6B 

18  Aug 

0650 

68.89 

18.3 

0.4 

7.5 

6C 

18  Aug 

0710 

68.89 

0.6 

Mb 

2.3 

7A 

18  Aug 

1101 

51.84 

24.0 

25.8 

26.8 

3  Units 

-2/3 

m 

:  for  the  multiplication  factor  are 

1  x  10"14. 

K-2  x  IQ"14 

2 

and  for  C„ 

N 

are 

Data  unavailable. 
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At  infrared  through  millimeter  wavelengths,  absorption  and  refraction 

2 

caused  by  water  vapor  become  the  dominant  contributors  to  .  Analytic 

applications  appropriate  for  infrared  and  mmw  propagation  are  credited  to  Hill 

et  al.  (1980)  and  Hill  and  Clifford  (1981).  Neglecting  secondary  effects  by 

other  atmospheric  constituents  (carbon  dioxide,  ozone,  aerosols),  they  derived 

expressions  for  the  log  amplitude  variance  of  infrared  through  mmw  radiation 

2  2 

propagation  as  functions  of  C.J,  ,  C  ,  and  C,^.  The  real  component  of  the 
refractive  index  for  a  propagating  spherical  wavefront  is  (Hill  et  al,  1980) 


c  2  _  at2  ct2  V2  Cl_ 

N  ■  T2  q2  Tq 


(2-13) 


where  A,j,  and  A^  are  frequency-dependent  dimensionless  sensitivity  coefficients 
that  account  for  changes  in  refractivity  and  absorption  due  to  temperature  and 
humidity  fluctuations.  McMillan  et  al  (1983)  define  AT  and  A^  for  mmw 
propagation  as 


A 


T 


77. 6P  1710q' 

- +  - 

It  t 


x  io  6 


(2-14) 


and 


A 


q 


1710q' 
.  T  , 


x  10  6 


(2-15) 


where  P  is  in  millibars,  q  is  in  grams  per  cubic  meter,  and  T  is  in  degrees 
Kelvin.  For  the  temperature-humidity  cross-structure,  Hill  et  al  (1988) 
suggest  that  r^q  is  unity.  Under  this  assumption,  Equation  (2-3)  shows  that 
C„  can  be  obtained  from 


(2-16) 
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A  detailed  justification  for  this  assumption  is  provided  in  Hill  et  al. 

2 

(1988).  The  Equation  (2-13)  components  for  mmw  CN  derived  from  DPG  data  are 

N  2  2 

presented  in  Table  9.  In  deriving  these  components,  C^.  and  were  obtained 
from  Tables  6  and  7,  and  was  calculated  using  Equation  (2-16). 


Comparison  of  the  mmw  components  in  Table  9  with  each  other  and  with 
the  optical  wavelength  provided  in  Table  8  illustrates  the  importance  of 

humidity  fluctuations  in  desert  terrain.  The  temperature  structure  term 

2  2 
(containing  C^.  )  in  Table  9  is  comparable  in  magnitude  to  the  optical  in 

Table  8,  which  is  due  almost  entirely  to  effects.  However,  the  humidity 

structure  and  humidity-temperature  cross-structure  terms  also  make  significant 

2 

contributions  to  mmw  .  The  sign  of  the  cross-structure  term  is  negative 
during  midday  trials  because  the  sign  of  is  negative  and  the  sign  of  is 
taken  from  r^,  which  is  positive  during  midday.  In  most  cases,  the  cross¬ 
structure  term  is  nearly  equal  in  magnitude  to  the  sum  of  the  temperature  and 
humidity  structure  terms.  Therefore,  if  rT  is  unity,  the  results  in  Table  9 
indicate  that  mmw  is  much  smaller  than  optical  during  the  day  and  much 
larger  than  optical  at  night.  For  comparison,  the  magnitudes  of  mmw  C^ 


'N 


with  r^  of  zero  and  unity  are  presented  in  Table  9.  The  diurnal  variation  of 

mmw  CN  appears  to  be  much  more  reasonable  under  the  assumption  of  r,^  equal 

to  zero.  As  noted  in  Subsection  2. A,  DPG  coherence  data  suggest  that  r_, 

2  1C* 

equal  to  zero  is  the  better  approximation  at  DPG.  Direct  mmw  CN  measurements 
were  not  available  to  verify  the  calculated  in  Table  9. 


2.6  COMPARISONS  OF  ACTINOMETER  AND  SCINTILLOMETER  DATA 


The  purposes  of  Phase  I  of  this  study,  the  demonstration  or  feasibility 

2 

phase,  were  to:  (a)  obtain  concurrent  radiation  and  CN  data  over  a  typical 
wintertime  snow  surface  and  (b)  develop  data  reduction  and  display  routines. 
Phase  I  was  originally  scheduled  to  begin  on  16  February  1988  and  last  for  2 
weeks.  One  week  before  the  planned  start,  DPG  had  a  7-  to  10-cm  snow  cover. 
With  temperatures  forecasted  to  be  above  freezing,  the  snow  cover  was  pre¬ 
dicted  to  dissipate  within  a  week.  In  order  to  make  measurements  over  snow, 
the  test  was  started  a  week  early  even  though  all  of  the  instrumentation  was 
not  operational.  Instrumentation  problems  occurred  throughout  Phase  I,  and 
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Table  9.  Components  of  Equation  (2-13)  for  Millimeter-Vave  CM  Computation. 
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Data  unavailable. 


the  test  vas  extended  to  4  weeks  in  order  to  collect  data  from  most  of  the 
instruments.  The  scintillometer  produced  valid  data  during  the  first  week  and 
part  of  the  second.  During  the  rest  of  Phase  I,  the  scintillometer  operated 
intermittently  because  of  a  series  of  problems  such  as  failure  of  the  LED  in 
the  transmitter. 

2 

An  analysis  of  relationships  between  available  Phase  I  scintillometer 
and  actinometer  radiation  data  was  performed  on  several  data  sets.  One  2-day 
data  set  was  collected  during  the  first  week  (10  and  11  February  1988). 
Temperatures  on  both  days  ranged  from  7  °C  during  the  day  to  -7  °C  at  night. 
Also,  skies  on  both  days  were  cloudy  early  in  the  day  and  then  decreased 
throughout  the  day.  Winds  were  strong  (5-20  m  s-^)  from  the  northwest  for 
much  of  10  February  and  light  and  variable  on  11  February.  Snow  cover  was 
10.2  cm  on  the  morning  of  10  February,  decreasing  to  2.5  cm  by  the  afternoon 
of  11  February. 

2 

Solar  radiation  and  measurements  for  10  and  11  February  1988  are 

shown  in  Figures  7  and  8,  respectively.  As  expected,  the  incident  short-wave 

radiation  on  10  February  rises  steadily  to  a  peak  around  solar  noon  and  then 

decreases.  The  few  dips  in  the  radiation  curve  represent  periods  of 

2 

cloudiness.  The  scintillometer  values  remain  steady  throughout  the  day 
and  then  start  to  decrease  about  3-4  hr  after  sunset.  readings  remain  low 

during  the  night  and  then  start  to  rise  again  about  1  hr  after  sunrise. 

2 

Little  relationship  was  found  between  C.t  and  incident  short-wave  radiation 

,  .  .  N  2 
during  the  daytime  on  10  February.  In  contrast,  on  11  February  the 

readings  tracked  with  changes  in  solar  radiation  (Figure  8).  There  are  two 

main  meteorological  differences  between  these  days:  (a)  the  10.2-cm  snow 

cover  of  10  February  had  almost  dissipated  by  11  February  and  (b)  the  wind 

speed  was  much  greater  on  10  February. 
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:ident  Short-Wave  Radiation 


An  interesting  Phase  I  case  occurred  on  16  February.  This  day  was  cloudy 

and  windy,  with  intermittent  snow  pellets  between  1200-1600  Mountain  Standard 

2  -12  -2/1 

Time  (MST).  Figure  9  shows  unusually  high  values  (10  m  )  between 

1200-1600  MST,  which  possibly  are  attributable  to  the  precipitation.  One 

2 

possible  explanation  for  the  high  measured  values  is  a  signal  reflected 

from  the  snow  pellets.  Another  possible  explanation  is  vibration  induced  on 

the  scintillometer  receiver  as  it  was  struck  by  the  snow  pellets.  Figure  9 

2 

indicates  that,  while  the  CN  values  peaked  during  precipitation,  the  incident 

N  2 
short-wave  radiation  curve  dipped.  When  the  precipitation  stopped, 

returned  to  normal  levels  and  the  radiation  curve  rose.  This  result  suggests 

that  further  testing  may  be  needed  to  determine  the  effects  of  various  types 

2 

and  intensities  of  precipitation  on  scintillometer  measurements. 


2 

Figure  9.  Scintillometer  CN  and  Incident  Short-Wave  Radiation  Data  for  16 
February  1988. 


Figure  9  also  illustrates  a  problem  that  began  at  approximately  1700  MST  on  16 
February  1988  when  the  scintillometer  began  to  give  a  constant  reading  of 
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10  m  .A  low  voltage  on  the  transmitter  was  found  when  the  instrument 
was  checked  the  following  day. 

The  Phase  I  objectives  of  demonstrating  the  measurement  capability  and 
developing  data  reduction  and  display  routines  were  partly  met.  All  of  the 
instruments  designated  for  the  Phase  I  test  were  deployed  and  data  were 
collected  from  them.  The  lessons  learned  during  Phase  I  were  helpful  in 
designing  the  August  Phase  II  trials.  Data  reduction  and  display  routines 

were  developed  which  made  data  interpretation  easier.  A  computer  program  was 

2 

written  to  convert  scintillometer  C,,  data  from  mV  to  the  commonly  used  units 
-2/3  N 

of  m  .  Another  program  was  developed  to  calculate  net  radiation  from  the 
four  actinometers .  A  computer  graphics  package  was  selected  to  produce  graphs 
for  analysis. 

2 

There  were  several  days  during  Phase  I  with  good  CN  and  incident  short¬ 
wave  radiation  measurements,  which  allowed  comparisons  to  be  made.  These 
measurements  did  not  correspond  well  when  there  was  snow  covering  the  ground. 
However,  there  was  only  one  day  during  Phase  I  with  the  desired  snow  cover. 

More  extensive  measurements  during  the  wintertime  with  snow  cover  are  needed  to 

2 

define  any  relationships  between  CN  and  other  meteorological  variables  during 

N  2 

that  time  of  the  year.  Better  correspondence  between  (?N  and  incident  short¬ 
wave  radiation  was  found  on  Phase  I  days  with  no  snow  cover.  Because  these 

2 

data  showed  promise  for  a  relationship  between  and  radiation,  this  possibi¬ 
lity  was  further  explored  during  Phase  II. 


The  scintillometer  was  operated  continuously  during  Phase  II  and  per¬ 
formed  well  with  no  problems.  Comparisons  between  Phase  II  scintillometer  and 

2 

actinometer  data  were  made  to  further  define  a  relationship  between  C.,  and 

N 

radiation.  Data  were  collected  from  0930  Mountain  Daylight  Time  (MDT)  on  16 
August  to  1200  MDT  on  18  August  1988.  Fifteen-minute  averages  were  used  to 
obtain  a  workable  data  set,  yet  not  smooth  out  fluctuations  due  to  transient 
changes  in  cloud  cover.  Meteorological  conditions  on  16  and  17  August  in¬ 


cluded  clear  skies  and  light  and  variable  winds  during  the  daylight  hours. 
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Temperatures  ranged  from  35  °C  at  midday  to  near  10  °C  at  night.  A  high,  thin 
cloud  cover  advected  over  the  test  site  after  sunrise  on  18  August. 


2 

Figure  10  shows  the  scintillometer  data  for  16-18  August  1988.  Peak 

values  of  C  ^  occur  during  the  afternoon  heating  as  extreme  surface  tempera- 
N  2 

tures  produce  large  density  gradients.  values  reach  minima  during  the 

transition  periods  when  the  atmosphere  becomes  quasi-adiabatic.  Although  the 
2 

magnitude  of  momentarily  drops  below  instrument  threshold  as  transition 

occurs,  the  data  from  this  event  are  averaged  with  the  data  obtained  before 

the  transition  to  produce  a  measurable  15-min  reading.  After  the  transition 

2 

minimum  at  sunset,  increases  again  during  the  nighttime  as  a  thermal 
imbalance  is  produced  by  radiational  cooling. 


TIME  (MDT) 

2 

Figure  10.  Scintillometer  Record  for  16  -  18  August  1988. 
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2 

It  should  be  noted  that  the  transition  period  minima  do  not  occur 

exactly  at  sunrise  and  sunset  because  of  the  thermal  lag  between  the  surface 

2 

and  air  temperatures.  The  minima  occur  15  to  30  min  after  sunrise  and  1 
hr  before  sunset.  These  results  are  similar  to  the  results  obtained  by 

2 

Walters  and  Kunkel  (1980)  at  White  Sands  Missile  Range  (WSMR),  where  the 
minima  were  about  1.5  hr  after  sunrise  and  0.5  hr  before  sunset.  Various 
factors  could  contribute  to  the  differences  between  DPG  and  WSMR,  including 
soil  type  and  height  of  measurement.  (The  measurement  heights  in  the  DPG  and 
WSMR  studies  were  2  and  8  m,  respectively.) 

Three  types  of  radiation  (incident  short-wave,  emitted  long-wave,  and 

2  " 
net)  were  compared  with  .  Incident  short-wave  radiation  versus  C^*"  for  16, 

17,  and  18  August  is  illustrated  in  Figure  11.  As  shown  by  the  figure,  the 
2 

variations  in  closely  track  changes  in  incident  short-wave  radiation 
during  the  daytime. 

2 

Emitted  long-wave  radiation  versus  is  shown  in  Figure  12.  The  best 

correspondence  between  the  two  is  found  during  the  daytime,  although  the 

2 

correlations  were  not  as  high  as  found  between  CN  and  incident  short-wave 

radiation.  Similar  analyses  of  the  transition  and  nighttime  data  indicated  no 

2 

interrelationships  between  CN  and  long-wave  radiation  during  transition 
periods  and  a  low  correlation  at  night. 

2 

The  third  radiation  parameter  compared  with  was  net  radiation.  Net 

radiation  was  computed  by  first  finding  the  net  short-wave  radiation  and  net 

long-wave  radiation  and  then  adding  the  two.  Net  short-wave  is  defined  as 

incident  short-wave  radiation  minus  reflected  short-wave  radiation,  while  net 

long-wave  radiation  is  defined  as  incident  long-wave  radiation  minus  emitted 

long-wave  radiation.  With  the  exception  of  the  empirically  derived  incident 

long-wave  radiation  (see  Subsection  2.2),  all  radiation  values  were  measured 

with  Epploy  actinometers .  As  shown  in  Figure  13,  the  correspondence  between 
2 

CN  and  net  radiation  was  again  high  during  the  daytime  and  not  significant 
during  the  transition  periods  and  at  night.  As  noted  earlier,  18  August  was 

the  only  day  during  Phase  II  with  some  cloud  cover.  It  is  interesting  to  note 

2 

the  correspondence  between  dips  in  both  the  CN  and  incident  short-wave  radia¬ 
tion  curves  on  this  day,  as  shown  in  Figure  14  between  0900  and  1200  MDT. 
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These  dips  are  periods  when  the  clouds  obscured  the  sun.  Similar  results  were 
noted  on  partly  cloudy  days  during  Phase  I. 

In  summary,  both  incident  short-wave  and  net  radiation  are  highly  corre- 
2 

lated  with  optical  CM  during  the  daytime  hours  of  clear  and  cloudy  days. 

N  2 
These  results  are  encouraging  for  our  goal  of  establishing  a  climatology 

using  available,  continuously  operated  instrumentation.  Unfortunately, 

results  were  not  so  good  during  the  transition  and  nighttime  phases,  and  it 

appears  that  another  method  will  have  to  be  found  for  these  parts  of  the 

diurnal  cycle.  No  Phase  II  testing  was  performed  during  periods  with  high 

winds  or  precipitation.  In  addition,  Phase  II  data  were  collected  over  only  a 

3-day  period.  A  longer  testing  period  is  needed  under  a  variety  of  meteorolo- 

2 

gical  conditions  in  order  to  verify  and  radiation  relationships. 
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Versus  Emitted  Long-Wave  Radiation  on  16  -  18  August  1988 
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Figure  13.  Optical  CN  Versus  Net  Radiation  on  16  -  18  August  1988 
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SECTION  3.  APPENDICES 

APPENDIX  A:  A  DESCRIPTION  OF  SCINTILLOMETER  OPERATION 


Scintillometer  operation  is  based  on  certain  assumptions  about  atmospher¬ 
ic  turbulence.  Consequently,  scintillometer  performance  is  dependent  on  how 
well  the  atmosphere  conforms  to  these  assumptions.  The  atmosphere  is  continu¬ 
ously  adjusting  to  the  flux  of  heat  by  turbulent  mixing,  which  creates  a 
spectrum  of  turbulent  elements  (eddies)  of  varying  densities.  Density  discon¬ 
tinuities  between  eddies  cause  changes  in  the  atmospheric  refractive  index.  A 
model  of  turbulence  initially  proposed  by  Kolmogorov  is  applied  to  the  refrac¬ 
tive  index  power  spectrum.  The  model  depicts  a  refractive-index  spectral 

2 

density  proportional  to  C.,  ,  representing  the  strength  of  refractive  index 

N  2 

fluctuations.  Tatarskii  (1961)  defines  C^  as  a  function  of  differences  in 
the  refractivity  (N)  at  position  x  and  at  scalar  separation  distances  (x+r) 
within  the  inertial  subrange. 


CN2  =  (N(x)  -  N(x+r))/r2/3  (A-l) 

where  the  overbar  indicates  a  time-averaged  quantity. 

The  Kolmogorov  model  is  assumed  to  represent  turbulence  in  the  inertial 
subrange,  a  range  of  eddy  sizes  from  a  few  millimeters  to  several  meters,  in¬ 
creasing  as  a  function  of  height  above  the  surface.  In  this  model,  turbulent 
energy  is  injected  at  the  larger  eddy  sizes  (the  outer  scale)  and  is  trans¬ 
ferred  to  smaller  scales  as  the  larger  eddies  break  up.  Energy  is  transferred 
to  progressively  smaller  scales  essentially  without  loss  until  the  dissipation 
range  (inner  scale)  is  reached.  In  the  dissipation  range,  turbulent  energy  is 
converted  to  heat  through  viscous  dissipation.  The  real  atmosphere  differs  in 
varying  degrees  from  the  ideal  turbulent  state  presented  in  this  model,  which 
assumes  the  existence  of  steady-state  uniform  turbulence.  Also,  density  dis¬ 
continuities  may  be  very  weak  in  near-neutral  conditions.  These  factors 
constrain  the  precision  and  range  of  measurements  made  by  wave  propagation 
techniques. 


A-l 


Scintillometers  are  designed  to  operate  using  turbulence  scales  within 
the  inertial  subrange.  The  principle  of  optical  scintillometer  operation  can 
be  illustrated  with  a  simple  eddy  model  presented  by  Clifford  et  al.  (1974). 
Figure  A-l  shows  an  optical  path  (OA),  established  between  the  transmitter  and 
receiver,  and  a  turbulent  eddy  of  radius  r  at  path  position  C-  The  turbulent 
eddy  represents  a  density  discontinuity  moving  across  the  optical  path,  and 
its  index  of  refraction  creates  a  new  path  OTA.  This  eddy  is  illuminated  by  a 
spherical  wave  of  wavelength  A.  The  illuminated  eddy  produces  a  diffraction 
pattern  observed  at  A,  the  size  and  contrast  of  which  are  determined  by  C,  r, 
A,  and  the  refractivity  fluctuation  AN  of  the  eddy.  Eddies  at  path  position  C 
which  satisfy  the  condition  OTA-  OA  =  A/2  produce  strong  scintillations  at  A 

through  destructive  interference.  This  condition  is  satisfied  by  a  Fresnel- 
zone  size  eddy  of  radius 

r  =  [ AC( 1-C/ A) ] 1/2  (A-2) 

Eddies  of  radius  larger  or  smaller  than  r  have  focal  lengths  shorter  or  longer 
than  A  and  therefore  contribute  less  to  the  variance  in  irradiance  received  at 
A.  Sensitivities  to  the  various  portions  of  the  optical  path  are  defined  by 
weighting  functions.  Weighting  functions  are  determined  by  appropriate 
transmitter-receiver  aperture  designs. 


Figure  A-l.  The  Geometry  of  a  Simple  Eddy  Model  Consisting  of  a  Source  0 

Illuminating  an  Eddy  of  Radius  r  at  Position  C  Along  the  Optical 
Path.  The  Resultant  Scintillations  are  Observed  at  Position  A 
(After  Clifford  et  al.,  1974). 
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2 

The  variance  in  received  irradiance  is  related  to  through  the 
Tatarskii  (1961)  first  order  theory  of  scintillation.  This  theory  is  based  on 
simple  superposition  of  the  effects  of  independent  eddy  encounters  along  the 
path.  It  predicts  that,  for  a  given  path  length  A  and  free-space  wavenumber  K 

(K=2n/A),  the  proportionality  relationship  between  the  log-amplitude  variance 

2  2 

in  irradiance  of  the  spherical  wave  (crt  )  and  is 


0.124  K7/6  A11/6 


(A-3) 


2 

Equation  (A-3)  is  valid  for  less  than  0.3.  In  stronger  turbulence,  the 

illuminating  wavefront  is  no  longer  spherical  due  to  distortion  by  previous 

2  2 

eddies,  and  the  proportionality  between  and  o  no  longer  remains 
constant.  This  condition,  known  as  saturation,  constrains  the  length  of  the 
optical  path  over  which  the  scintillometer  can  be  operated.  Saturation  and 
the  signal-to-noise  ratio,  which  decreases  with  path  length,  are  the  two  major 
constraints  on  instrument  range. 


A-3 


INTENTIONALLY  BLANK 
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APPENDIX  B:  DESCRIPTION  OF  DPG  SPECTRUM  ANALYSIS  PROGRAM 


Spectrum  analysis  of  time  series  data  provides  statistical  measures  of 
instrument  performance  and  insight  into  the  physics  governing  the  interactions 
between  sets  of  time  series  data.  Simple  linear  correlation  coefficients  are 
not  applicable  to  time  series  data  because  significance  tests  are  based  on  the 
hypothesis  that  the  variate  samples  are  uncorrelated.  Because  significant 
autocorrelation  often  exists  in  time  series  data,  this  hypothesis  is  seriously 
violated.  Also,  the  true  relationship  between  two  time  series  can  be  obscured 
by  an  in-phase  relationship  at  some  frequencies  and  an  out-of-phase  or  phase- 
lagged  relationship  at  other  frequencies.  A  solution  to  some  time  series 
analysis  problems  lies  in  the  multivariate  analyses  applied  to  spectra  of  time 
series.  These  analyses  can  be  accomplished  after  a  Fourier  transformation  of 
the  series  from  the  time  domain  to  the  frequency  domain.  Transformation  par¬ 
titions  the  power  (energy  expenditure  per  unit  time)  of  the  time  series  into 
harmonic  frequency  components.  For  a  given  variable,  "power"  is  equivalent  to 
variance  in  that  variable  per  unit  time.  The  total  power  of  the  process  is 
equal  to  the  sum  of  the  contributions  by  the  harmonic  components.  The  advan¬ 
tage  of  power  spectrum  analysis  over  the  usual  multivariate  analysis  is  that 
the  amount  of  power  contributed  by  one  harmonic  is  independent  of  the  ampli¬ 
tudes,  phases,  and  frequencies  of  the  other  harmonics  in  the  time  series 
(Koopmans,  1974). 

Relationships  between  two  time  series  can  be  evaluated  by  examination  of 
the  spectra  from  the  two  series  and  their  complex  products,  the  cross  spec¬ 
trum,  phase  (PHASE),  and  coherence  (COH).  The  cross  spectrum  is  represented 
by  the  cospectrum  (COS)  and  quadrature  (QUAD).  The  cospectrum  is  analogous  to 
an  in-phase  covariance  between  the  two  spectra,  and  the  quadrature  spectrum  is 
a  similar  measure  phase-shifted  1/4  wavelength  (90°).  COS  and  QUAD  define  the 
cross  spectrum  covariance  components  in  a  Cartesian-like  coordinate  system, 
while  PHASE  and  COH  express  these  components  in  a  normalized  polar  coordinate 
system.  PHASE,  a  representation  of  the  angular  relationship  between  two 
spectra,  is  defined  as 


PHASE  =  ATAN2( QUAD/COS), 


(  B— 1 ) 


where  ATAN2  is  defined  in  FORTRAN  as  the  arctangent  function  expanded  to  the 
range  -n  to  rt.  PHASE  is  presented  in  degrees.  COH,  the  squared  coefficient 
of  coherence,  is  a  measure  of  the  correlation  between  two  time  series  as  a 
function  of  frequency  given  by 

COH  =  (COS2  +  QUAD2)  /( SPCTR1 )*( SPCTR2 ) ,  (B-2) 

where  SPCTR1  and  SPCTR2  are  the  spectral  components  of  time  series  1  and  2, 
respectively.  COH  is  dimensionless  and  ranges  in  magnitude  from  0  to  1. 
Statistical  significance  tests  can  be  applied  to  coherence  data  as  detailed  in 
Koopmans  (1974). 
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Fourier  transformation  in  the  DPG  spectrum  analysis  program  is  achieved 
using  a  modification  of  the  Fast  Fourier  Transform  (FFT)  procedures  provided 
to  DPG  by  the  NOAA  ERL/VPL,  as  described  in  Kaimal  and  Gaynor  (1983).  The 
spectra  produced  by  this  FFT  are  logarithmically  smoothed  and  scaled  to 
meter-kilogram-second  (MRS)  units  by  multiplying  each  harmonic  component  by 
its  frequency.  Table  B.l  is  a  sample  of  spectra  and  their  derived  components. 
Spectral  analysis  techniques  are  discussed  in  detail  by  Koopmans  (1974)  and 
Jenkins  and  Vatts  (1968). 


Table  B.l.  Spectra  From  Trial  1A,  Channel  5  (T)  and  6(q)  Data. 

SPECTRA  DATA  SUMMARY 


SPECTRAL  VALUES 


Frequency 

HZ 

Spectr 
(T'T' ) 

Spectr3 

(q'q' ) 

COS 
(T'q' ) 

QUAD 

(T'q') 

PHASE 

Deg 

C0H 

0.0020 

0.0175667 

0.0000034 

-0.0001853 

-0.0000027 

181. 

0.582 

0.0029 

0.0144672 

0.0000011 

-0.0001148 

-0.0000127 

186. 

0.825 

0.0039 

0.0869024 

0.0000079 

-0.0007883 

0.0001364 

170. 

0.927 

0.0049 

0.0414242 

0.0000036 

-0.0003323 

-0.0000642 

191. 

0.766 

0.0063 

0.0456963 

0.0000045 

-0.0004011 

0.0001053 

165. 

0.832 

0.0083 

0.1778915 

0.0000133 

-0.0015203 

-0.0000497 

182. 

0.974 

0.0107 

0.0733853 

0.0000055 

-0.0006052 

0.0000972 

171. 

0.926 

0.0141 

0.1611553 

0.0000137 

-0.0014603 

0.0001042 

176. 

0.968 

0.0189 

0.0975598 

0.0000063 

-0.0007142 

0.0000711 

174. 

0.844 

0.0248 

0.0953864 

0.0000075 

-0.0007909 

0.0000843 

174. 

0.887 

0.0320 

0.1434356 

0.0000110 

-0.0011914 

0.0000648 

177. 

0.903 

0.0417 

0.0889597 

0.0000073 

-0.0007702 

-0.0000329 

182. 

0.916 

0.0543 

0.1153469 

0.0000086 

-0.0009583 

0.0000238 

179. 

0.923 

0.0708 

0.1138459 

0.0000079 

-0.0009086 

0.0000170 

179. 

0.914 

0.0920 

0.1053497 

0.0000070 

-0.0008173 

-0.0000068 

180. 

0.908 

0.1197 

0.0912712 

0.0000064 

-0.0007125 

0.0000264 

178. 

0.873 

0.1550 

0.0768547 

0.0000060 

-0.0006317 

0.0000307 

177. 

0.867 

0.2010 

0.0639370 

0.0000048 

-0.0005009 

0.0000179 

178. 

0.822 

0.2606 

0.0552272 

0.0000044 

-0.0004310 

0.0000193 

177. 

0.769 

0.3380 

0.0502979 

0.0000038 

-0.0003805 

0.0000116 

178. 

0.764 

0.4387 

0.0428086 

0.0000035 

-0.0003198 

0.0000213 

176. 

0.687 
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APPENDIX  C:  LIST  OF  SYMBOLS 


A  =  refractive  index  coefficient  for  optical  wavelength  (78.7  x  10~^  K  mb-^) 
Aq  =  dimensionless  mmv  sensitivity  coefficient  for  humidity,  ^1710q ^-6 

A_  =  dimensionless  mmw  sensitivity  coefficient  for  temperature, 

♦  ^3)x  ID'6 

B  =  Bowen  ratio  (dimensionless)  CpP  w'T'/(eLv  w'e') 
c  =  calories 

2  -32  -2/3 

0^  =  humidity  structure  parameter  ((g  m  )  m  ) 

-3  -2/3 

^Tq  =  temperature-humidity  cross-structure  parameter  ((K-g  m  )m  ) 

2  -2/3 

CN  =  refractive  index  structure  parameter  (m  ) 

Cp  =  specific  heat  at  constant  pressure  (0.24  c/g/K) 

2  2  -2/3 

CT  =  temperature  structure  parameter  (Km  ) 

°C  =  degrees  Celsius 

e  =  water  vapor  pressure  (mb) 

eg  =  saturation  vapor  pressure  (mb) 

_2 

g  =  acceleration  due  to  gravity  (9.8  ms  ) 

H  =  surface  layer  sensible  heat  flux  (K  m  s  *) 

He  =  latent  heat  flux  (mb  m  s-^) 

Hs  =  H  derived  from  scintillation  techniques  (K  m  s “*) 
k  =  Von  Karman  constant  (0.4,  dimensionless) 

K  =  absolute  temperature,  Kelvin 

- - 

L  =  Obukhov  length,  meters  (-u^  T/(kg  w'T')) 

Ly  =  latent  heat  of  vaporization  (586  c  g~*  at  20  °C) 
m  =  meter,  or  prefix  milli- 
n  =  real  part  of  the  refractive  index 
N  =  refractivity  (dimensionless,  N=(n-l)xl0^) 
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P  =  pressure  (mb) 

q  = 

%  : 

R  = 


r  = 

rTq 

s  = 

T  = 
T 


* 

t  = 

u  = 

V  = 
w  = 

V  = 

u*  = 

2  = 

2  = 
O 

a  = 

t  = 

V  = 

e  = 

K  = 

X  = 

A  = 

P  = 

2 


_3 

absolute  humidity  (g  m  ) 

_3 

:  humidity  scaling  factor  (w'q'/u^,  g  m  ) 

-2 

net  radiation  (W  m  ) 
scalar  separation  distance  (m) 

2 

=  temperature-humidity  cross-structure  correlation  coefficient  (C^,  /C^C^) 

second 

temperature  (K) 

:  temperature  scaling  factor  (w'T'/u^,  K) 
time  (s) 

horizontal  wind  speed  (m  s  *) 
volt 

vertical  wind  speed  (m  s 
watts 

friction  velocity  (m  s 
height  above  the  surface  (m) 
roughness  length  (m) 

dimensionless  wavelength-dependent  structure  parameter  coefficient 

ratio  of  molecular  weights  of  water  vapor  and  dry  air  (0.622) 

water  vapor  correction  factor  (dimensionless) 

potential  temperature  (K) 

free  space  wavenumber  (1/m) 

wavelength  (m) 

optical  path  length  (m) 

-3 


air  density  (g  m  ) 

=  log  amplitude  variance  in  irradiance 


-1. 


=  square  root  of  the  variance  in  vertical  wind  (ms  ) 
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